Hydrostatic equilibrium equation and Newtonian limit of the singular f(R) gravity 
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I. Introduction 



A. J. Bustelo and D.E. Barraco 
FaM.A.F., Universidad Nacional de Cordoba 
Ciudad Universitaria, Cordoba 5000, Argentin^\ 

We derive the hydrostatic equilibrium equation of a spherical star for any gravitational Lagrangian 
density of the form L = y / —gf(R). The Palatini variational principle for the Helmholtz Lagrangian 
in the Einstein gauge is used to obtain the field equations in this gauge. The equilibrium hydrostatic 
equation is obtained and is used to study the Newtonian limit for f(R) = R — jg. The same 

procedure is carried out for the more generally case f(R) = R — -^p^ a fi t giving a good Newtonian 
limit. 

PACS numbers: 04.50.+h,98.65.-r,98.70.Vc 
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On From recent studies it seems well established that our universe is currently in an accelerating phase. The evidence 
of cosmic acceleration has arisen not only from the high-redshift surveys of type la supernovae [l], 0] , but also from 
the anisotropy power spectrum of the cosmic microwave background [1, 3 ■ One of the most accepted explanations 
is that the universe has been dominated by some form of dark energy for a long time. However, none of the existing 
^vq [ dark energy models are completely satisfactory. 

■ It is possible to find other explanations for cosmic expansion using field equations other than Einstein's equations. 
0^ Recently, some authors have proposed to add a i?" 1 term in the Emstem-Hilbert action to modify general relativity 
Q. They have obtained the field equations using the second order formalism, varying only the metric field, and have 
thus derived the so-called fourth-order field equations. Although the models were obtained using corrections of the 
Hilbert-Einstein Lagrangian of type R n , where n can take a positive or negative value to explain both the inflation at 
an early time and the expansion at the present timera. Apart from this ad hoc justification, there are also theoretical 
motivations for non- linear gravity from M-theory Q|. However, this model still suffer violent instabilities \A which 
maybe deleted by the addition of a i? 2 -like term to the Lagranngian [§J or by account of quantum effects Q. A big 
review on various modified gravities and its applications to cosmology was analyzed in and the Newtonian limit 
of these fourth order theories has been studied by Dick 

There are also theories that are obtained from a Lagrangian density L t = y/—gf(r) + L m , which depends on the 
scalar of curvature and a matter Lagrangian does not depends on the connection, on this Lagrangian we can apply 
Palatini's method to obtain the field equations [111, [Hj]. In references [U, EH, we showed the universality of the 
Einstein equation using a cosmological constant. More recently, M. Ferraris, M.Francaviglia and I. Volovich published 
the same result [HI, [Tj] . For these theories we have studied the conserved quantities [l5| , the spherically symmetric 
solutions [l6| , the Newtonian limit [Til HH, [l7| , and the Cosmology described by FRW metric [l8| . 

In [l7| . it was shown that is very difficult to test these models in the (post-) Newtonian approximation. The reason 
for this is that the departures from Newtonian behavior are both very small and are masked by other effects, due to 
the fact that these departures have to be measured when the body is moving "through" a matter-filled region. 

Recently, Vollick [l9[ used the Lagrangian f(R) = R— f^, together with the Palatini variational principle to prove 
that the solution of field equations approach the De Sitter universe at a late time. This result was obtained using, 
the well known property of the vacuum solutions of these theories, that in a vacuum, or in the case of T = constant, 
the solutions are the same as in General Relativity with a cosmological constant, even when f(R) is not analytical at 

fl = op3,[IE[l2|. 

On the other hand, solutions corresponding to different cosmological constants are allowed by some of these theories. 
Therefore, the homogeneous and isotropic vacuum solution for these theories is the deSitter space-time with different 
cosmological constants, except when one of the allowed cosmological constants is A = 0, which corresponds to flat 
space-time. Thus, the inclusion of 1/R curvature terms in the gravitational Lagrangian provides us with an alternative 
explanation for cosmological acceleration. The generalization to the case of scalar tensorial theories was analyzed in 
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In the present work, we use the first order or Palatini formalism [llj in order to obtain the hydrostatic equilibrium 
equation, for theories that are obtained from the Lagrangian density L = yJ-gf(R) which depends on the scalar of 
curvature, and a matter Lagrangian that does not depend on the connection. This is done by means of a Legendre 
transformation which in classical mechanics replaces the Lagrangian of a mechanical system with the Hclmholtz 
Lagrangian [2l|. From this equation we have found the Newtonian limit in the particular case f(R) = R — #g, is a 
good Newtonian limit in same cases and is not in others. This limit has been treated in the literature [H, Hjjwith 
no consensus about what the background metric is, that has to be perturbed. In this sense, the difference from 
previous works is that we obtain this limit from the hydrostatic equilibrium equation and we assume nothing about 
the background metric. In addition, we generalize the above result for the case f(R) = R - -^p 1 9 ^-. Also, for this 
theory we see that the departure from Newtonian behavior is smaller as n increase. 

In the next section we review the field equations and the basic structure of the theory. In section III we derive 
the static spherically symmetric solution. In section IV we derive the equilibrium hydrostatic equation, modeling the 
star as a perfect fluid for f(R) theories of gravity. In the following section we determine the weak field limit from 
this equation. In section VI we calculate the weak field limit for singular f(R) gravity model, and in section VIII we 
show the conclusions. This work is part of a graduation thesis [HJ 



II. Helmholtz Lagrangian in the first order formalism 

The action for an f(R) gravity in the Jordan(original) gauge with the metric g ab is given by [2f 

' ' ^f(R)d A x + S m (g ab ^), (1) 



2a r , 

Here y/—gf(R), is a Lagrangian density that depends on the curvature scalar R = R ab (T de )g ah , S m is the action for 
matter represented symbolically by ip and independent of the connection, and a m — 8nG, where we have taken c = 1. 

We consider that we have a smooth one-parameter (A) family of field configurations starting from given fields g ab , 
T° de and ip, with appropriate boundary conditions and, denote by Sg ab , ST ab , Sip the corresponding variations, i.e., 

Sg ab — ^j^\\ = oSX, etc. Variation of the action Sj with respect to g a b yields the field equations [ll[ 

f'(R)R ab - \f{R)g ab = a m T ab , (2) 

where f'(R) = (Ms) an d the dynamical energy-moment tensor of matter is generated taking variation of the action 
of matter with respect to the Jordan metric tensor: 



SS n 

The trace of field equation @, gives 



\ J T ab ^d 4 xSg ab . (3) 



f'{R)R-2f(R) = a m T, (4) 

where Q fixes a unique relation, in general non-linear, between R and T. Varying the action with respect to the 
connection, and recalling that this is fixed al the boundary of hi , we obtain: 



K(.V=gg aD f(R)) = o, (5) 

this equation define, the connection F which is not the metric connection. The metric connection satisfies the equation 
V c g a b = 0. From the field equation it follows that the connection T is given by the Christoffel symbols of the 
conformal transformed metric: 

g ab = f'(R)g ab . (6) 

The metric g ab defines the Einstein gauge in which the connection T is metric compatible. Therefore, instead of using 
we will use V c . 

The action fl]) is dynamically equivalent to the Hclmholtz action, given by [2l[ 

S H = — I d 4 x^[f^(p))+p(R(g) - 0(p))] + 5 m (. 9ab », (7) 
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where p is a new scalar variable, and <j)(P) is an auxiliary scalar field determined by the following equation 

0/, 



OR 



\r=<P(p) = p. (8) 



Taking variations of the action with respect to p, and assuming that the Lagrangian of matter does not depend on 
the connection, we have: 

0(p) = R(g ab ). (9) 
From the form of the scalar of curvature in the Jordan gauge and from (|6|) we obtain 

<f> = Rf{4>) = Rp, (io) 

where R = R a b(J^ C d e )g ab is the scalar of curvature of the metric g ab . 

Transforming (J7J) to the Einstein gauge this action becomes the standard Hilbert-Einstein action with an additional 
scalar field 

2a m J p p z 



Se = — I d 4 x^-g[R-^ + + S m (p- L g ab ^)}, (11) 



and choosing 4> as the scalar variable (note that 4>(p) = R ab (T)g ab i.e. is the scalar of curvature in the Jordan gauge) 
leads to 

Se = ^~ [ d i x^—g[R-2V( ( f)]+S m {g ab [f{4>)]-\i>) (12) 



2a 

where V(4>) is the potential: 

m = w • (13) 

Now the variation of the action with respect to g ab gives the field equations of the action ([1]) in the Einstein gauge 

Rab - -^Rgab = V{(t>)g ab . (14) 

Then, considering IjlOp and taking trace of the before field equation we obtain 

- 2/(0) = a m T([/'(0)]~ 1 g afc )/'(0), (15) 
with the last equation being equivalent to equation ([¥]) but written in the Einstein gauge. 

III. Static spherically symmetric solution 

For the two cases of f(R) theories, those which emerge from second order formalism and those which emerge from 
first order formalism, the static spherical symmetric solution for vacuum has been studied in 

EE HI 

Let us work 

in the Einstein gauge and consider now a static spherically symmetric spacetime manifold with metric g ab in the 
standard form: 

d~s 2 = e air) dr 2 + r 2 duo 2 - e^dt 2 , (16) 

where the coordinates are curvature coordinates of the metric g ab . 

The nonzero components of the Einstein tensor [27| . G ab — R ab — ^Rg ab , are 

e~ a 1 

n6 _ f«p _ - af 7,rr , 7,r , 7,r <V <V7,r , 

^e-^ v -e (_ + T + _-_-^), 
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«_e-°(l-ra, r ) 1 



G\ = V - ^ (17) 

Let us now regard T t * and as assigned then the solution of the field equations (fT4|) . in the region r < ro where 
ro is the proper radius of the star, are 



(18) 



'"'/'"'Si 13 ''" 1 ' 1 * 1, (19) 



where M(r) is given by: 



MM - 1 - ff 2 { ^El _ MHi )dP (20 ) 

The spherical symmetry of the metric is, of course, invariant under the conformal change ©.Thus, to return to the 
original gauge, we have to perform the inverse transformation ^ and we can write the metric g ab in the form 

f {r) f (r) f'{r) 

It is convenient to put the metric in the standard form in order to know the physical meaning of the coordinates, 
therefore we change the coordinates from the curvature coordinates of the g a b to the curvature coordinates of the 
original metricg a b. In this coordinates r' = , r and the original metric g ab is 



ds 2 = A(r')dr' 2 + r' 2 dw 2 - B(r')dt 2 , (22) 



where 



l y/Fw dr ' 



Air') = = , (23) 

K J 1 2M(r>) ' V ; 



and 



M(r') is: 



flM = m [1 - ^P- + eM I fe °m {r - - T!)d ' % m 



iVTWJo V'W ' TO)] 2 

This results are similar to that obtained in [1(1 , in that work was also proved that in all cases the exterior metrics 
match correctly with the interior solutions. 



IV. The hydrostatic equilibrium equation 

In the original spacetime, Jordan gauge, the energy momentum tensor corresponds to a perfect fluid model: 

Tab(gcd) = pu a u b + p(g ab + u a u b ). (26) 
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In order to obtain T a b(f g c d) we note that the energy momentum tensor will be of the form of a perfect fluid 
because u a u a = u a u a = —1; i.e. u a Ub — f'(<p)u a Ub- Thus u a = —J — u a and 

T a b{g c d[f ifi)]' 1 ) = pu a Ub + p{g a b + u a u b ), (27) 

where p = p/ f '{(/>) and p = pj /'(</>)■ In order to derive the equation of hydrostatic equilibrium the field equation that 
we are interesting in, is 

fir n , \ 1 »m? "mP <£/'(</>) " / W / oeA 

Gr = — (1+r7 ' r) -^ = 7w" m = [7W" W ' (28) 

and making algebraic manipulations to separate 7 iT . we obtain 

_ 2M(r)+r 3 [ am p[/(0)]- 1 -y(0)] 
r(r-2M(r)) 



where M(r) is given by (JSOJ) and 1/(0) by (fT3| . 

The use of a conformal transformation , Einstein gauge, is only used as a trick which allows us to simplify calcula- 
tions, namely to calculate 7 jr . 

The matter action must be invariant under diffeomorphisms and the matter fields satisfy the matter field equations, 
thus [13: 

V b T ah = 0. (30) 

On the other hand, in the Einstein gauge the energy-momentum tensor is not covariantly conserved [281 ] - Therefore, in 
the Jordan gauge, a test particle will follow the geodesies of the metric connection but this is not true in the Einstein 
gauge. 

We assume the Jordan gauge as the physical, in order to research observational consequence of this hypothesis. 

The equation (|30p is written in Jordan gauge, and the energy momentum tensor is given by (|26[) since we model 
the star as a perfect fluid. Choosing the four-velocity pointing into the future like direction and normalizing it to 
u a u a = —1, it becomes 

u M = (vC4p), 0,0,0), (31) 

For our metric, in Jordan gauge (|2ip . the only nontrivial component gives, 

dp _ (p + p) tt dg tt , , 

dr> ~ 2 9 dr> ' {6Z) 
and in order to combine this equation with equation (|29p . we now have to consider equation ([6]), then: 

dp ( P + P )d 1 (p + P ) <y f{4>) 



dr' 2 dr' ^/fWl dr ' 

Making a change by taking into account r' = , = 



(33) 



and putting everything in terms of r' on the right hand side of equation (|29[) . we can write: 



gr =fl . t d^fW) )( r'H am pirm- i -rwvm + 2M(r / ) 

dr' v / 7 7 P) dr ' r'(r'-2M(r')) 

Finally, we are ready to write the equilibrium hydrostatic equation in the Jordan gauge: 



1 ifflw 

/~F77~7T Hi > X 



Cp + p)dr' 2^ ^T^y d r ' 



( r /3 (a m p[/ / (0)]" 1 -r(0)^(0)) + 2Af(rO 1 d^fW) m 

r'{r' - 2M{r')) dr ' 

The equation (|36[) can be seen as equivalent to Tolman-Oppenheimer-Volkoff equation, (for a different approach see 
[29|). It is easy to show that for General Relativity, i.e. f(4>) = (f>, the equation (|3li>)) becomes the known result. This 
equation is very difficult to solve, for any model different from General Relativity, even if we consider an incompressible 
fluid. 
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V. Newtonian limit of the f(R) gravity 

In order to obtain the Newtonian limit of the f(R) gravity from the hydrostatic equilibrium equation, we assume 
p <C p, Gm(r) <C r, where m(r) = Air p(x)x 2 dx, and take p small enough so that it produces a good approximation 
in the following Taylor expansion 

/(0(T)) = f(T) « /(0) - p^ | T=0 +0{p\ (37) 

taking into account T = 3p — p ~ —p. Similar expansions were made for the other quantities appearing in (|36|) . In 
particular an expansion of V(<p) is given by the equation (I13p 



vm~-WL- 0( ^i±mmi^ mf"m\ T =o , 



Thus, taking into account equation ([JJ in order to solve -^\t=o, and considering this equation in the case T = we 
can write: 



Now, let us approximate the integral equation (|25p . First we consider its integrand: 



^.wya.^^tort. (40) 



/'(*) ' 2[/'(*)] 2 [/'(0)P 2[/'((0)] 

The integral part of M(r') can be written in the form 



7'(0 4 2[/w ; 7 M/'(o)F 6 [/'(op 

where we have assumed p continuum. Finally, the expression (|25[) takes the form 

( } T(^r ( + 7W)^ T=a) + v7W) ' ( } 

where m(r') is the usual integral expression for the mass. Another term in the expression (|36[) is 

aw»[/ (0)] - / « [ W) (1 + J7(o) ^It=o) + ^|t=o] - 5 (43) 

The terms in the brackets in (|43p can be ignored because has terms of superior order. Plugging the approximations 
into equation (|36|) and working consistently with the conditions in this limit gives 



dr'~/'(0) 1 r' 2 j 6/'(0) r d r ' ' 1 > 

solving the derivative in the right hand side of equation (|44[) and making one more approximation keeping to linear 
order in p we have: 

_ * ~ g _ / + P /"(0) gm ^ f45) 

dr /~//( ) r* 6/'(0) r 2 /'(0) (0(O)/"(O)-/'(O))dr' - 1 j 

If the particular choice of theory of gravity is such that /(0) = 0, we have a theory without cosmological constant, 
and the last equation can be written in the form: 

dp ^ p Gm(r') pa m /"(0) dp 

dr' ~ /'(0) r' 2 2 [f'{Q)Y { dr' h 1 ' 
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The last expressions give the weak field limit from the hydrostatic equilibrium equation for the theory, and they are 
totaly in accordance with previous works [23j. However, the Taylor expansion (|37[) can be done in general, when p is 

2 

small enough, but not for singular theories of type f(R) = R — f^, in which the condition in order to keep the linear 
term in p is pa m -C a [22j], where a is essentially the vacuum energy density with p vac = a/a m ~ 10~ 31 [gr/cm 3 }. It 
is obvious that the mean density of our solar system p ~ 10~ n [gr/cm 3 ] does not satisfy this condition, and then we 
can not use the above approximation. Therefore, in the next section we will study the weak field limit for this kind 
of theories in the case of p being of order of the mean density of our solar system. 

VI. Newtonian limit of the singular f(R) gravity 

Theories like those presented in [lj|, i.e. R— 2 ^, which have a pole in R = 0, seem to explain the present accelerated 
expansion. But, there are still some obscure points about how to obtain their behavior in a weak field limit, particularly 
what metric we have to use as a background metric. Dominguez and Barraco p3| perturbed around the de Sitter 
space-time, which is the maximal symmetric vacuum solution of the theory. Sotiriou [22j claimed that this procedure 
was erroneous, when p is near the mean density of the solar system, which, according to the final part of the last 
section, is correct. Therefore, he perturbed around the flat metric, which is not a solution of the theories but it could 
be near the vacuum solutions. 

In our present work we suppose nothing about the background metric in order to obtain the Newtonian limit, and 
this is one of the principal difference with the previous works. First we find this limit for the theory: 

m = R ~M- (4?) 

For this theory, equation Q gives 



<f>(T) = -(-a m T ± V« 2 m T2 + 4a2) , (48) 

where we have to choose the minus sign in front the square root, in order to return to General Relativity when a — > 0. 
We consider, p < p, r' » Gm(r') and: 

— « 1. (49) 

OL m p 

Here we are considering [22] p — 10 -11 [gr /cm 3 ], a ~ 10 _67 [ey] 2 so | |~ 10~ 21 where T w —p. Then we can 

state: 

4>(T) k, a m p + (—)a, 
a m p 

2 a 

f{4>) w a m p+ -( )a, 

6 a m p 

3 a m p 

f"(<f>)~-l( — ) 2 — - (50) 
6 a m p a m p 

Next, we have to approximate the others quantities appearing in the hydrostatic equilibrium equation, and take 
into account that in this approximation M(r') « Gm(r'), where m(r') is the ordinary mass appearing in General 
Relativity. Then, plugging into (I36p we have, 

dp Gm(r') 1. a ,,dT_ Gm(r') . 

_ T7 = P 71^ - o ) T7 1 r 1 ) 

ar r 6 a m p dr r 
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Gm(r') 1 a 2 d{3p - p) 
^^^"st^ — A3 - " (51) 

we have already neglected the pressure but we can't neglect its derivative. Therefore, the second line is obtained from 
the first by substituting T — 3p — p. Using (|4"5)) we have, in the weak field approximation: 



P _, 2 +?( — ) 317 . ( 52 ) 



dp Gm(r') ^ 1 , a ^ 2 
dr' r' 2 3 a m p dr' 

where the terms of order a 2 are kept only in gradients of p, where they are relevant. The extra term in equation (|52[) 
is the correction of lowest order in to the purely Newtonian behavior. Therefore, the last equation gives a good 
Newtonian limit for the theory. This limit is in accordance with the result obtained by Sotiriou. 
We can show more generally that the theories 

1 a n+1 

f(R) = R- * , (53) 

n + 2 R n 

also has a good Newtonian limit. Equation ([J| for these theories give 

<t> n+1 + a m T<f> n - a n+1 = 0, (54) 

where the last equation has an approximate solution of the form <p(T) rs —a m T « a m p. Then using the expressions 

2 a " 

f{4>) ~a m p+ — — ( )a, 

n + 2 a m p 

n + 2 a m p 
n + 2 a TO p a m /9 

considering these approximations in the context of equation (|36|) and carrying out exactly the same procedure as 
before: 

_ dp_ _ Gm(r') n (n + 1) a , n+1 jp , , 

dr'~ P r' 2 2{n + 2ya m p' dr r [ ' 

Where the last term is the correction of lowest order in -^-^ . Note that for the case n = 1, the last equation returns 
to equation (f52"|) and for n = 0we return to General Relativity 

Finally, the above equation shows that these theories have also a good Newtonian limit. 

VII. Conclusions 

We have obtained the hydrostatic equilibrium equation in the Jordan gauge, where the matter obeys simple conser- 
vation laws. This equation is the generalization of the Tolman-Oppcnhcimcr-Volkoff equation in General Relativity. 

We have studied the Newtonian Limit, in this gauge, from the hydrostatic equilibrium equation, and in section V 
we found an expression in the weak field limit for the pressure gradient. In this section we saw that this result is not 
correct if we consider f(R) = R — gravity models with density typically related to a Newtonian regime. However 
the expansion (|46|) for -j4 is perfectly accurate for this singular gravity if the conditions are such that a m p <§; a. And, 
if this is the case, the theory does not give a good Newtonian Limit as previously shown Barraco and Dominguez 23] . 

In section VI we have shown that the singular theory, f(R) = R— and its generalization, f(R) = R — ^^^tt-j 
have good Newtonian limits for densities such as the mean density of the solar system. The difference with previous 
works [221.I23T] is that we found the Newtonian limit, in the spherically symmetric case, with no assumption about the 
background metric. Finally, in the general case, we have shown that n controls the departure from purely Newtonian 
behavior. 
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